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ABSTRACT 

The  aerodynamic  characteristics  of  oblique  wings  in  an  inviscid, 
incompressible  flow,  linearized  for  small  wing  camber  and  incidence, 
is  studied  under  the  assumption  that  the  wing  aspect  ratio  Al,  is  high. 

Apart  from  the  addition  of  a dominant  upwash  correction  of  the  order 
/H,  resulting  from  the  sweep  of  the  center  line,  the  present 
analysis  differs  from  the  classical  lifting-line  theory  in  that  the  flow 
field  next  to  the  wing  section  (the  inner  solution)  is  affected  by  a 
component  of  the  wake  vorticity  parallel  to  the  center  line,  and,  hence, 
is  not  locally  two-dimensional.  A crucial  aspect  of  the  analysis  involves  the 
behavior  of  the  three-dimensional  correct  ions  near  the  leading  and  trailing 
edges,  which  require  special  attention,  lest  nonuniformities  arise.  The 
results  determined  from  matching  the  inner  and  outer  solutions  exhibit  a 
strong  asymmetrical  spanwise  influence  of  the  wake  vorticities,  with  a 
lift  increase  on  the  downstream  wing  panel  and  a lift  reduction  on  the 
upstream  panel.  Results  obtained  are  compared  with  surface-lift  distri- 
butions generated  by  an  inversed  method  for  yawed  elliptic  planforms, 
and  with  span  loadings  generated  by  a panel  method  for  elliptic  flat 
plates  (wings  with  zero  camber)  as  well  as  an  ESP  (extended-span  planform) 
wing.  For  in  the  range  of  10  to  20,  good  agreement  in  the  comparison 
is  consistently  found,  and  the  improvement  over  the  strip  (local  2-D)  theory 
is  shown  to  be  great.  Recast  into  a rational  fraction  (in  a form  similar 
to  that  used  originally  by  Prandtl),  results  obtained  can  be  improved 
further  and  shown  to  be  adequate  for  aspect  ratio  down  to  at  least 
4.33  corresponding  to  a 5:1  ellipse  at  30°  yaw.  The  report  also 
furnishes  computed  (finite-part  of  the)upwash  data  which  will  be  useful 
in  other  related  subsonic  and  transonic  applications. 


1.  INTRODUCTION 


The  effect  of  far-wake  vorticities  is  essential  in  any  analysis  of 
three-dimensional  ( 3 “D ) lifting  surfaces.  In  Prandtl 's  (1918)  lifting- 
line theory  of  a high-aspect-ratio  wing,  the  calculation  of  this  effect 
is  greatly  simplified  by  replacing  the  wing  with  a line  of  singularity 
comprising  all  bound  vortices.  A restriction  in  Prandtl 's  work,  as  well 
as  the  more  systematic,  asymptotic  theory  of  Van  Dyke  (196*0  and 
Ashely  and  Landahl  (1965),  is  the  assumption  that  the  wing  center  line 
(a  reference  curve  in  the  theory)  is  straight  and  unswept  (being  perpen- 
dicular to  the  main  flow).  The  present  study  considers  extensions  of 
this  approach  to  planar  wing  probl ems* i n which  the  afore  mentioned 
restriction  is  removed.  The  primary  purpose  of  this  note  is  to  point  out 
certain  distinct  features  of  such  a theory,  concentrating  mainly  on  the 
development  for  a steady  incompressible  potential  flow  past  an  oblique  wing 
(of  which  the  center  line  may  be  taken  as  being  straight).  Comparison 
with  results  of  more  exact  calculations  by  other  methods  will  be  made. 

The  work  reported  represents  a special  case  of  the  theory  of  high- 
aspect-ratio  wings  involving  curved  center  lines;  while  the  analysis  is 
simplified  owing  to  the  absence  of  the  center-line  curvature,  a few  of 
the  basic  features  brought  out  below  remain  crucial  ingredients  in  the 
more  general  theory  to  follow.  An  abr  idged  vers  ion  of  this  work  has  been 
made  in  the  form  of  a technical  note  (Cheng  1978).  Apart  from  a fuller 
presentation  and  discussion  of  the  theory,  this  report  gives  additional 
comparisons  of  the  theory  with  the  panel  method  in  the  lower  aspect-ratio 
range  and  for  an  ESP  oblique  wing  not  previously  considered. 
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Among  earlier  methods  employing  the  lifting-line  idea  to  swept  wings, 
(see,  for  example  , Jones  and  Cohen,  1957),  the  most  well  known  is,  perhaps, 


that  of  Weisinger  (19**2),  which  fails,  however,  in  the  limit  of  an 
infinite  aspect  ratio.  Solutions  to  elliptic  lifting  surfaces  at  yaw 
have  been  given  early  by  Krienes  (19^0),  based  on  a superposition  of 
products  involving  Lame's  functions.  Readers  are  referred  to  Jones  and 
Cohen  (1957)  for  a helpful  delineation  of  Krienes*  analysis,  where 
applications  of  certain  results  for  unyawed  elliptic  planforms  are  made. 
Krienes  treated  the  direct  problem  for  an  el 1 i pt i c flat  plate,  and 
obtained  results  for  an  5:1  ellipse  at  0-30°  yaw.  Of  the  five-term 
truncated  series  used  therein,  three  terms  were  symmetric  spanwise,  it  is  not 
clear  if  the  remaining  terms  would  adequately  delineate  the  asymmetric 
span  load  of  interest.  Generalization  of  the  lifting-line  theory  for 
a wing  with  a curved  center  line,  as  well  as  a wing  in  side  slip,  has  been 
made  by  Dorodnitsyn  (19AA),  who  noted  the  significance  of  the  logarithmically 
large  upwash  due  to  yaw.  However,  the  analysis  involved  an  ad_  hoc  treatment 
in  which  the  lifting  surface  is  replaced  by  a lifting  line  at  the  quarter- 
chord  location  and  the  resulting  upwash  correction  is  evaluated  at  the 
three-quarter-chord  location.  The  result  was,  in  any  case,  restricted  to 
small  departures  from  a straight,  unswept  center  line. 

A case  of  curved  center  line  has  been  studied  by  Thurber  (1965)  who 
considered  a crescent-moon-shaped  wing  from  the  view  point  of  asymptotic 
methods;  but  the  inner  solution  and  the  matching  problem  were  not  considered 
(the  upwash  calculation  also  contained  errors).  Oscillating  high-aspect- 
ratio  wings  with  curved  center  lines  have  been  treated  by  the  author 
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(Cheng  1976);  however,  upwash  calculation  for  the  low-frequency  and 
quasi-steady  cases  has  not  been  given,  and  the  formal  inner  solution 
presented  therein  is  incomplete  (see  below).  Apart  from  the  hope  of 
gaining  physical  insights  and  simplicity  for  analyzing  transonic  oblique 
wings  and  animal  swimming  propulsion  (Jones  1972,  1977;  Lighthi 1 1 
1975;  Chopra  and  Kambe  1977;  Cheng  1976),  the  present  asymptotic  approach 
has  been  motivated  by  a desire  to  implement  the  current  computer-oriented 
3-D  methods  which,  though  very  powerful,  are  by  no  means  free  from  both 
cost  and  storage  limitations  (see  Bauer,  Garabedian,  Korn  and  Jameson 
1974,  Jameson  and  Caughey  1977,  Ashley  and  Rodden  1972).  Problems  involving 
curved  center  lines  in  steady  and  unsteady  flows  will  be  treated  fully  in 
separate  papers.  Extension  to  transonic  flow  problems  involving  nonlinear 
component  flow  have  been  studied  by  Cheng  and  Meng  (1978,  1978b). 
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2.  PRELIMINARY  REMARKS 


In  Prandtl's  theory,^  the  3~D  effects  appear  mainly  as  a local 
incidence  correction  and  can  be  calculated  readily  as  one  half  of  the 
upwash  at  the  wing  trace  in  the  Trefftz  plane.  This  simplicity  is  lost, 
if  the  center  1 i ne  has  a nonvanishing  sweep  angle  or  curvature;  the 
solution  in  this  case  cannot  be  unambiguouslydetermined  without  recourse 
to  a more  systematic  analysis  employing  matched  asymptotic  expansions. 

Certain  features  of  the  theory  absent  in  the  classical  works  are  noteworthy 
and  may  be  inferred  from  rather  elementary  considerations.  As  will  be  seen 
below,  most  of  these  features  may  be  traced  to  the  presence  of  a spanwise 
component  of  the  local  wake  vorticity. 

Logarithmic  upwash 

The  most  obvious  among  these  features  is  perhaps  the  dominant 
induced  upwash  resulting  from  the  sweep.  From  the  view  point  of  an  asymptotic 
theory,  this  additional  upwash  is  significant  in  that  it  adds  terms  of  the 
order  fit~'2n  fit,  to  corrections  which  are  otherwise  of  the  order  fit., ' , 
where  fit,  is  an  appropriately  defined  aspect  ratio  to  be  given  later.  As 
pointed  out  by  Cheng  (1976),  the  wake  vorticity  (using  Prandtl 's 

notation)  shedded  local ly  has  a component  along  the  swept  center  line, 
which  gives  rise  to  a logarithmically  singular  upwash  near  the  center 
(lifting)  line.  Therefore,  this  raises  the  magnitude  of  the  upwash  in  the 
inner-flow  region  next  to  the  wing  by  a factor  of  . 

The  near-wake  influence 

With  a nonvanishing  spanwise  component  of  the  wake  vorticities,  the 
flow  region  next  to  the  wing  section  can  no  longer  be  well  represented  by 
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wakeless  2-D  divergence-free  flow  as  in  most  c 1 ass  ica  I ana  lyses . Thus,  in 
addition  to  the  upwash  induced  by  the  far  wake  (to  be  determined  by 
matching),  the  inner  solution  must  properly  account  for  the  upwash  at  the 
upper  and  lower  airfoil  surfaces  induced  by  this  component  of  the 
vorticity  right  behind  the  trailing  edge  (which,  to  the  leading  order, 
is  proportional  to  the  local  value  of  s i nA  -d  ) . 

Problems  of  edge  singularities:  nonuniformities 

As  in  most  classical  lifting-surface  theory  dealing  with  a direct 
problem,  in  which  a finite  upwash  (slope  of  the  camber  surface)  is  prescribed 

over  the  wing,  it  is  required  that  the  singularity  of  the  velocity  at 

-•/*  y. 

the  leading  edge  and  trailing  edge  be  no  stronger  than  p and  p , 
respectively.  In  above,  p stands  for  distance  from  the  edge.  For 
wing  upwash  which  becomes  logarithmically  infinite  at  the  edges  like 
p , it  will  be  required  that  the  nearby  velocity  field  behaves  no 
worse  than  l*i  p . At  the  leading  edge,  these  requirements  may  be  alterna- 
tively replaced  by  one  requiring  i ntegrabi 1 i ty  (with  respect  to  distance) 
of  the  surface  pressure  or  surface  speed.  At  the  trailing  edge,  the 
requirement  is  equivalent  to  the  Kutta-Joukowski i condition. 

In  constructing  analytic  solutions  to  the  reduced  problems  in  the 
asymptotic  theory  for  large  , it  is  essential  to  adhere  to  the  foregoing 
conditions  at  successive  orders  of  /R(,the  requirements  tantamount  to 
enforcing  un i form  val id i ty  of  the  expansion  at  both  edges.  This  is  so 
because  partial  derivatives  with  respect  to  distance  along  the  center 

line,  say  a/*/s  appear  in  certain  equations,  which  would  raise  the 
order  of  a.i  algebraic  or  logarithmic  singularity,  unless  the  location 
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of  the  singularity  is  independent  of  y'  (see  Section  *4  below).  Homogeneous 
solutions  with  proper  singular  behaviors  must  be  added  when  necessary 
to  fulfill  the  stated  requirements  and  to  avoid  nonuniformity. 


We  note  in  passing  that  an  additional  logarithmic  upwash  appears 
if  the  relative  wing  motion  is  unsteady,  or  if  the  center  line  curvature  is 
comparable  to  the  reciprocal  of  the  span. 


3.  COORDINATES  AND  DESCRIPTION  OF  THE  FULL  PROBLEM 

The  following  analysis,  employs  two  right-handed  Cartesian 
system.  In  (x,y,z)  system,  the  reference  wing  plane  is  z = 0 and 
makes  an  angle  of  attack  oL0  with  the  undisturbed  flow;  thus,  the 

velocity  vector  at  the  infinity  appears  in  this  system  as  (U  cos  , 

0 

0 , I*  sin  ),  where  U is  the  free-stream  speed.  The  second  system 
(x'.y'.z1)  is  generated  from  the  first  by  a rotation  about  the  z-axis 
(cf.  Fig.  1). 

x'z  coiAx- Sin  Ay  , y'-SinAK+cosAy,  z = z'  (3.1) 

The  y'-axis  is  taken  as  the  center  line  of  the  oblique  wing.  Note  that 

the  wing  planes  lies  in  Z/  = * = 0 and  that  A.  is  the  sweep  angle.  The 

variables  x*  and  z'  are  normalized  by  the  root  chord  CQ  (measured  parallel  to 

x'axis);  y',  along  with  x,y  and  z,  are  normalized  by  the  half  span  b. 

The  primed  system  is  employed  to  describe  the  inner  region  near  the  wing 

section  and  the  unprimed  one  for  the  outer  region. 

The  perturbation  potential  <f>  of  the  full  linearized  problem  must 

satisfy  the  3~D  Laplace  equation,  being  regular  everywhere  except  in 

approaching  the  wing  and  the  trailing  vortex  sheet;  it  must  also  vanish 

at  the  infinity  except  near  the  far  wake.  In  approaching  the  upper  and 

lower  surfaces,  z = z (x,y),  the  linearized  impermeability  condition, 

w 


transferred  to  the  wing  plane  z = o,  is 

(it)  = (jlL*  . * 


(3.2a) 
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and  in  approaching  the  trailing  vortex  (TV)  sheet,  the  continuity  of 
pressure  and  normal  velocity  requires 


ttt£Dr7  ttfUL- 0 


(3.2b) 


where  the  subscript  w and  TV  signify  the  wing  and  the  trai 1 ing-vortex 
sheet,  respectively,  and  the  double  bracketjj^  j^stands  for  differences 
across  the  wing  or  the  TV  sheet.  The  velocity  is  allowed  to  be  unbounded 
at  the  leading  edge  but  must  be  integrable;  the  pressure  or  ^ on  the 
wing  is  required  to  remain  bounded  at  the  trailing  edge  so  that  the 
Kutta  condition  is  satisfied. 

We  observe  in  passing  that,  the  free-stream  velocity  vector  could 
be  alternately  set  to  be  (U,  0,  0),  along  with  the  omission  of  — 
from  Eq.  (3.2a)  (note  that in  this  alternative  system  differs  from 
that  in  the  original  (x,y,z)  system  by  -at#).  This  alternative  system 
with  the  free-stream  velocity  as  (U,  0,  0)  has  been  used  in  Cheng  (1978) 
and  is  equivalent  to  the  (x,y,z  ) system  adopted  in  the  present  work,  for 
all  practical  purposes.  The  present  (x,y,z  ) system  with  (U  cos  oL9  , 0, 

U sinct*  ) is  preferred^as  it  avoids  a conceptual  difficulty  which  one 
would  encounter  in  the  transfer  of  the  wing  boundary  condition  to  the 
wing  plane  (z  = 0)  in  the  case  of  a high-aspect-ratio  wing  when  the 
angle  of  attack  is  comparable  to 

Throughout  the  development  in  Sections  k and  5,  an  aspect  ratio 
defined  as  ft,  = * *> A.  is  used.  Note  that  the  distance  from  tip  to  tip 
2 £ is  less  than  the  span  2b  = 2 l • cos  A . Another  aspect  ratio 
*„  = i l/c9  = A,  secA  will  be  used  later  for  the  convenience  in 
comparing  results  for  the  same  wing  at  different  sweep  angles.  The  aspect 


ratio  A,  is  preferred,  as  the  theory  developed  requires 


'/lb  = « 1 


which  is  significantly  different  from  Ca/2  l = /R~'  « /.  The  parameter 
controlling  wing  camber  and  incidence,  including  the  angle  oia  , is  denoted 


by  fc  . The  wing  ordinate  normalized  by  £ C0  , is  written  as 
Z(x' ,y') . 


EPftr, 


*4.  INNER  PROBLEM 


The  linearized  problem  under  a small  £ admits  a development  for 
high  aspect  ra^io  (for  finite  x'.y'.z') 

♦' 2 = *• r *'*■  * + " (',  l> 

allowing  a weak  logarithmic  dependence  of  $ and  <Pt  on  . The  first 
two  terms  satisfy  the  2-D  Laplace  equation  in  x‘  and  z*.  The  conditions 
on  the  wing  and  the  vortex  sheet,  Eqs.  (3-2),  through  the  transform  Eq.  (3.1), 


yield 


(£)=£.*  -*>  |l fro. 


va  *n 

-yP  = o i (A. 2 a , b ,c) 


(H)  - Zmk  2 • It^I  .I’srl = 0 (4-3 


where  m si  tan  A , and  o i0  = - °^/e;  2m[|^fJJ  results  from  the  component 

of  trailing  vorticity  parallel  to  the  center  line. 

The  solution  ^ satisfying  the  boundary  conditions  (*4.2),  fulfilling 

the  edge  conditions  mentioned,  and  yielding  a zero  disturbance  velocity  at 
2 2 

large  x*  + z1  , is  provided  by  the  classical  2-D  thin  airfoil  theory.  As 
it  is  well  known,  the  voriticity  strength  []>*./»*-]]  on  the  wing  is  a 
solution  to  the  integral  equation  ( 

J-  py.  ( .jfrftfoj.  - 2-z-Z 

27T  9k'  ° 


x;  -x' 


where  P.V.  signifies  the  Cauchy  principal  value;  a1  and  b'  locate  the 
leading  and  trailing  edges,  respectively. 
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The  solution  <^)  which  meets  the  same  edge  singularity  requirements 
and  allows  matching  with  the  outer  solution  can  be  obtained  as  the  sum 

Pt  = f y>p  +-  vT*'  (**-5) 

where  V^°  is  an  anticipated  upwash  correction  to  be  determined  by 

n P 

matching  later;  y is  a (particular)  solution  satisfying  the  (nonhomogeneous) 
TV  boundary  conditions  in  Eq.  (A. 3b),  and  f is  a (wakeless)  2-D  thin 
airfoil  solution  with  its  wing  upwash  so  chosen  to  make  pt  fulfilling 
Eq.  (4.3a).  The  solution  *pP  can  be  obtained  as 

= - ZmR.p.  ($  ?W/dr)  t %P  (4.6) 

where  R.P.  signifies  the  real  part,  £ is  the  complex  variable  x‘  + i z ' , 

W is  the  complex  potential  pg  + i ifo  ; is  a wakeless  harmonic  function 

(not  considered  in  Cheng  1978)  , needed  to  keep  hence  <£»  , from 

becoming  more  singular  than  <pg  at  the  edges  (see  below).  To  fulfill 
Eq.  (4.3a),  the  jump  in  across  the  wing  must  satisfy  the  same 

integral  equation  as  Eq.  (4.4),  with  the  R.H.S.  replaced  however,  by 

- f Zm  x'^Z - ('7Vh/'9x')v  t 

+ £ ( Z - ZlB  - *2.  a')  - vr.  <“-7> 

In  the  absence  of  the  Kutta  condition,  the  , together  with  the  part 

of  p contributed  by  -(  vVS  /sz'  )w  in  Vt  above  would  belong  to  the 
family  of  e i gensol ut ions  for  <pf  . 

In  choosing  p P , one  notes  that  the  singular  behavior  of  (p  is 
dominated  by  the  first  member  of  pP  in  Eq . (4.6)  as  2ma  1 ( 9 /“dy' 
at  the  leading  edge  and  2mb  ' (' * at  the  trail ing 

edge.  For  cases  wherein  the  camber  slope  is  nowhere  infinite,  it 
suffices  to  take 


12 


{pf=  Z”a'lp,  * Z"Jc'E  ^ R p.  [i  f 


(A. 8) 


where  E is  the  coefficient  in  ^ Ax'~  E(b'-x')*.  For  the  case  with 
infinite  camber  slopes  at  the  edges,  the  prescription  is 

tPHP  = 2ma'^<p0Le  t-  2.»rb'ff,<t>o£  ^-9) 

where  is  a solution  of  the  ^ type,  with  its  x'-derivative  tending 
to  at  L.E.  and  vanishing  no  less  rapidly  than  (b'-x1)  at  T.E.; 

is  similar,  with  supercript  L.E.  replace  by  T.E.  and  a1  by  b‘.  An 
important  example  for  this  case  is  logarithmically  infinite  cambers 
supporting  finite  pressure  jumps  at  the  edges. 

In  terms  of  the  pressure  coefficient  based  on  the  dynamic 
pressure  of  the  component  flow,  the  pressure  jump  across  the  wing  can  be 
obtained  from 


[M  = ~lr  A.  " h-*  " ZmX  9/ix’^o  f9y% 


The  span  loading  is  PU  times  the  potential  jump  at  the  trailing  edge  at 
the  y-(spanwise)  station  and  can  be  computed  from  the  jumps  in  , <f 

and  if?  at  the  trailing  edge  at  the  corresponding  y'-station;  in  normalized 

H 

form 


s.l.s  * I<>IL = 


r i 


b' 


v 

<#>'a/  (nr)'1  jTCy)  C tan'(*,/e')  + 7 Sj*  z'J  -(**)  f Xcp»l J *•  d*>  R-P-(lG  ) * 
t XC'v'm  iy*  - *»’«.  *■[*-"(»■/-) 'fspfj* 

t V,"  *'  » O ( ) , 


(4.12) 


where  |T  stands  for  , and  ff/')  = IT//'.)  ♦ ff  f + %P  I 


5.  OUTER  PROBLEM  AND  MATCHING 


* 


In  flow  region  removed  from  the  wing  sections,  the  spatial  variation 
is  scaled  by  b » ca  . and  the  solution  may  be  represented  to  the  leading 
order  by  concentrated  vortices  with  circulation  7^(y)  bounded  to  the 
center  line,  together  with  the  trailing  vortex  sheet.  The  perturbation 
potential  for  this  outer  problem  is  then  i - f,  * • • • • 

with 


(*'/■*>  = 


j_  r'  cw 

4rr  l, 


+ 


x-  my, 

iu-*y,)x4  (y-ydx+  z 


(5.1) 


anticipating  the  determination  of  7J  by  matching,  both  and  ^ in 
Eq.  (5-1)  are  normalized  by  BC0  1/  (same  as  for  ff1  and  f ' ) . 

In  approaching  the  center  line  (£  = y-»»y  — +0 , 1 -*o)  the  integrand 
shown  in  Eq.  (5.1)  becomes  singular  at  u = - / = 0 like 


gc  ) u ) 


rw * rt,y  * r , , (-•>«  i <5.2) 

u 1 + z * *u  -* 


where  Ru  = J ( K 64  + i u + f . By  substraction  of  terms 
shown  above  from  the  integrand,  the  resulting  integral  then  has  a 
(finite)  limit  as  £ and  t,  vanish.  This,  together  with  the  quadrature 
of  g,  gives  the  inner  limit  sought 


$ ^ tin  '("  2 ) + S$»  z J - 


(5.3) 


J.  * ( To  (y,)  - TcCy)  r.  m ( , .T  , 

+ » i - —y  L ' - v +•  * 

(5.3) 

+ o (cs+t)2-,  s'*:1)  t 

where  Oi  $ M,  ) is  expected  to  arise  from  Ai,  £ , . In  terms  of  the  inner 
variables  x'=  z'\ceSA  $ , ^'  = £‘/R,  * , and  observing  that  y)  = 

H(  , Eq.  (5-3)  can  be  expressed  as 

$ jjp  *7^  (y'c*SA)‘[  ^'(jr)  * 7 s9”  z ']  ~ 

-i»  *>'*'  iy  £*»*'(•?)*  1*9**']  * (5.4a) 

+ w fc/  ap  Me/  + 

♦ **"'*'  r.(o)ZCy)  * 0(  *,'V,  z",  <9 

where  % (y)  is  the  upwash  function  in  y (not  /') 

Z^5  7^'Otf  ^ WO  * 

- «i  ^ l> 1 fr  I - *'•"«  (l 4/- 1 is- 1 )]  * 

****  [,  — ^yl-y  ^ ' [ 1 " Jji.  rf/>  (s.'tb) 

and  £(/;=  7»(y) / Vo  = *T,  /<ty  • 

.1 

Eqs.  (4.12)  and  (5.4)  are  valid  expansions  of  in  the  common  domain 
/ « If  I « *.  , where  the  inner  and  outer  solutions  can  be  matched, 
provided  one  makes  the  identifications 

it  o'csa)  =ict/o , V“=  tnwzc/).  <s-5) 

The  upwash  correction,  hence  the  inner  solution  to  the  order  /ff  * , is  now 
determined . 

It  is  useful  to  point  out  that  is  the  finite  part  of  the  local 

induced  flow  angle  ot.  divided  by  £ , and  that2I  is  /Kol-  divided  by  the 

4 14 

sectional  lift  coefficient  at  y = 0.  The  first  term  in  Eq.  (5.4b)  cor- 
responds to  the  logarithmic  upwash  mentioned  earlier  which  dominates  the 
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finite  part  of  the  induced  velocity,  increasing  the  upwash  on  the 

aft  panel  (where  dVo/^O)  and  reducing  it  on  the  forward  panel  (where 

dTl/jy  > 0 ).  For  f^-distri  but  ions  of  interest,  dJl/dy  behaves  near 
-Vt 

y =±1  like  ( I-/2)  , this  together  with  the  logarithm  involving  4?, 

and  (I-/*),  leads  to  a 2 maximum  near  the  downstream  tip  at  I-  y=  O(K') 
where  the  induced  upwash  arises  from  0(^,  i* ) to  0(  A,  ).  There  is 

also  a minimum  in  2 near  the  upstream  tip  at  1 + y = 0(j4?(  ),  where  the 
induced  downwash  has  a magnitude  of  0(  At/*  ) . Although  the  magnitude  of 
the  induced  velocity  becomes  infinity  at  the  tips,  2 reverses  its  sign 
and  vanishes  on  each  wing  panel  at  a span  station  extremely  close  to  the  tips. 
This  tends  to  provide  a reasonable  description  of  the  span  loading  near 
the  tip,  in  spite  of  the  local  breakdown. 


6.  ILLUSTRATION  OF  THE  UPWASH  RESULT 


Figure  2 illustrates  the  spanwi se  £ d i str i but  ion  for  an  elliptic 
load  J*m  =y  I - / 1 and  an  ESP  load 

f.  = jrp  - (6.D 

The  latter  distribution  has  been  adopted  in  an  oblique-wing  design  study 
by  Black,  Beamish  and  Alexander  (1975),  and,  pertains  to  the  optimum 
span  load  for  minimum  induced  drag  with  given  lift  under  a fixed  wing- 
root  bending  moment  (Jones  1950).  For  a bending  moment  less  than  that 
of  the  elliptic  load,  the  wing  must  have  a longer  span.  A planform  based 
on  this  load  (assuming  a zero  camber)  has  been  referred  to  as  an  extend- 
span  planform  (ESP).  The  ESP  load  shown  in  Eq.  (6.1)  and  referred  to 
in  Fig.  2 corresponds  to  a particular  ratio  of  the  extended  span  to  the 
original  span  taken  as  1.15.  Results  are  shown  for  the  swept  angle 
A = 0 , 22.5°  and  **5°  with  a particular  aspect  ratio 

Zl/Ct  = 8.390$.  (6.2) 

S i nee  c©SA  appears  only  through  1*^  in  > cf.  Ecl-  (S.^b),  the 

graphs  can  be  used  for  other  aspect  ratio  by  simply  adding  to  2L 

-Fi  (6.3) 


* 

It  is  useful  to  note  that  if  a length  scale  C#  =^o/c  other  than  the  root 
chord  c0  is  used  for  the  inner  solution,  no  alteration  is  needed  in  the 
above  analysis,  except  changing  /R,  to  <T/^,  • The  relation  £ = **■>;/ C*~y-  0 

is  unaffected. 
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Figure  2 shows  that  the  departure  from  the  unyawed  upwash  distribution 
is  always  anti symmetr i c (in  y ),  and  is  considerable  larger  in  the  elliptic 
cases  (broken  line).  Interestingly,  the  degree  of  spanwise  nonuniformity 
in  / is  much  greater  for  the  ESP  load  (full  line),  varying  from  0.4  at 
the  center  to  1.10  at  y = 0.9.  This  is  because  the  unyawed  result  of  Jj 
for  the  ESP  load  is  itself  highly  nonumiform  spanwise.  In  passing,  it  may 
be  noted  that  the  exact  % for  ESP  for  A=  ° is  a linear  upwash  distribution 
(Jones  1950) 

Z « e,i  c,. |,| 

(marked  in  thin  dashes, in  Fig.  2),  which  is  recovered  reasonably  well  by 
a numerical  procedure  applying  to  the  ESP  load  of  Eq.  (5.6).  The  discrep- 
ancy shown  indicates  the  level  of  accuracy  to  be  expected  in  the  subsequent 
spanload  calculation  for  the  ESP  case  for  A.%0. 
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7.  COMPARISONS  WITH  OTHER  SOLUTIONS 


In  the  following,  comparisons  are  made  of  the  present  analysis  with 
exact  solutions  derived  from  an  inverse  method  and  with  results  from  a 
panel  method. 

7.1  Uniformly  Loaded  Oblique  Wings 
Wing  upwash  (slope)  supporting  a given  lift  distribution  can  be 
determined  with  relative  ease;  for  a uniform  distribution,  the  task 
reduces  to  evaluating  a line  integral.  Upwash  calculation  has  been  made 
for  several  planforms  with  uniform  load  in  this  manner  to  furnish  a basis 
for  assessing  the  theory.  Using  the  upwash  data  as  an  input , D is 
computed  by  the  present  theory  and  then  compared  to  the  exact  (uniform) 
load.  This  provides  a test  of  the  theory  as  a dj_rej:t  method.  Of  all 
examples  considered  with  Alm=ll/c0  = 10  to  AO  and  A = 0 to  1»5°,  a 
useful  feature  is  that,  over  most  Darts  of  the  wing,  the  upwash  can  be 
very  closely  represented  by  (with  c'  = b1  - a') 

V0(*’r)=  (A<r)+ B(;0~)  CV)  (7.D 

where  B and  C are  generally  small  and  decrease  with  increasing  aspect 
ratio.  This  fact  is  utilized  to  expedite  the  computation. 

Figure  3 presents  chordwise  lift  distributions  at  several  spanwi se 
stations  computed  from  such  an  upwash  over  a 20:1  elliptic  planform  at 
1*5°  yaw,  for  which  the  exact  Hcp'1  has  a normlized  value  unity  (a  similar 
result  has  been  presented  in  Cheng's  (1978)  note  for  a different  yaw  angle, 

A = 22.6°).  Good  agreement  with  the  exact  value  (unity)  is  seen  at 
the  inboard  stations.  This  represents  a great  improvement  over  QCp  *])  from  the 
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strip  theory (lower  graph)which  has  an  error  of  typically  5 to  20%. 

Near  the  tips  ( y / t>  -*■  £ 1),  the  lifting-line  solutions  deteriorates, 
as  expected.  At  the  80%  station  (dash  curves),  signs  of  breakdown 
appear  near  the  L.E.,  which  is  partly  caused  by  the  inadequacy  of  curve 
fitting  based  on  Eq.  (17)  near  the  tip.  However,  the  span  loading  (not 
shown)  appears  to  remain  satisfactory  even  at  the  80%  span  station. 

Other  examples  with  a lower  aspect  ratio  ( 10),  as  well  as  non- 

uniform  lift  distributions,  have  been  studied;  except  for  a slightly 
larger  discrepancy  near  the  tips,  similar  conclusions  can  be  drawn  for 
those  cases.  We  shall  return  to  the  question  on  applicability  of  the 
present  work  to  wings  of  lower  aspect  ratioslater  in  Section  7. 3- 

7.2  Oblique  Flat  Wings 

Vortex- 1 att i ce  and  wing-panel  methods  have  proven  adequate  in  many 
linear  lifting-surface  problems, t*  but  their  application  to  oblique  wings 


does  not  appear  to  exist  in  the  literature.  The  data  by  pane  1 method  used 
in  the  following  comparison (provided  by  Mr.  Ronald  Smith)  are  generated 
from  a computer  program  at  NASA  Ames  Research  Center,  based  on  an  exten- 
sion of  Woodward's  (1973)  method  applied  originally  to  symmetric  planforms. 


Figure  *4  presents  a comparison  in  span  loading  for  an  elliptic  flat-plate 
wing  of  A0  = 16.78  with  a straight  ^0%  chord  line  at  A5°  yaw.  For 
elliptic  flat  plates,  the  theory  gives  a span  loading  in  terms  of  elementaary 


functions 

S.L.  /a  cos  A 


= /T7‘  - * SmASin'y)  * 

+ o-il)  - 


- C$CA  [bl  i+3in/i I - 0-SihA)ln/cosAj]J 


(7.2) 


See,  for  example,  Ashley  and  Rodden  (1972)  and  Woodward  (1973). 


r 
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where  oL  is  the  angle  of  attach  and  k is  the  straight-axis  location  in 
fraction  of  a chord.  As  clearly  shown,  the  difference  between  the 
present  theory  (in  full  line)  and  the  panel  method  (in  open  circles) 
is  small  at  most  stations  as  compared  to  either  of  their  differences 
from  the  (uncorrected)  strip  theory  (in  dash).  The  numerical  data  from 
the  panel  method  are  generated  from  a run  employing  100  panels  over  20 
span  stations.  The  program  is  a part  of  the  NASA  Ames  wing-body  program 
for  a linearized  compressible  flow;  in  the  computation  the  free-stream 
Mach  number  was  set  equal  to  M^  = 0.10,  with  an  expected  compress i 1 i ty 
effect  comparable  to  one  percent. 

A similar  comparison  is  presented  in  Figure  5 for  an  ESP  flat  plate 
mentioned  in  Section  6,  with  same  A o and  A as  in  the  last  figure. 

Except  at  one  station  near  mid-span,  where  the  panel  method  gives  a 
slightly  higher  value,  the  agreement  is  seen  to  be  even  better.  (According 
to  a second  set  of  data  from  the  panel  method  for  the  same  problem  based 
on  10  span  stations  with  10  panels  for  each  span  station,  the  noticeable 
discrepancy  near  the  mid-span  disapears). 

7.3  Applications  Involving  a Lower  Mi 

The  usefulness  of  an  asymptotic  theory  based  on  » 1,  such  as 
the  present  one,  will  depend  on  its  adequacy  in  applications  where  4^, 
is  not  so  large.  As  it  is  quite  well  known,  Prandtl's  success  in  coping 
with  wings  of  lower  aspect  ratios  lies  largely  in  a special  solution  form 
retained  in  his  original  (i9'®)  work.  In  the  case  of  an  elliptic  flat 
plate  at  zero  yaw,  for  example,  the  lift  coefficient  is  obtained  and  commuted 
according  to  the  familiar  form  as 

ct=  znd./[,  + 2/^  ] (7.3) 
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where  M.  s(¥/rr)  A^a  , As  pointed  out  by  Van  Dyke  (196*0.  the  above  result 
is  consistent  with  the  (second-order)  asymptotic  theory  for  high  /R.  t 

CL~  ZTT<*  [l  ~ ] , (7. A) 

but  is  clearly  superior  than  the  latter  in  the  low  /R.  range.  Note  that, 
as  yR-*0,  Prandtl's  formula  in  rational  fraction^Eq.  (7-3)>yields  a limit 
which  is  larger  than  the  correct  slender-wing  limit  by  a factor  of  2,  but 
Eq.  (7.*0  gives  an  infinite  lift  coefficient  which  is  certainly  much 
worse.  For  ^-5-10,  or  , differences  of  Eq.  (7-**)  from  (7.3),  and 

from  the  more  exact  calculation^are  considerable.  The  accuracy  of  the 
present  theory  is  expected  to  be  no  better  than  its  counterpart  at  zero 
yaw,  Eq.  (7-*0»  and  some  very  definite  discrepancy  may  therefore  appear 
for  aspect  ratio  below  ten. 

Figure  6 presents  an  example  in  the  lower-aspect-ratio  range  in  which 
results  of  span  loadings  of  a 5J1  elliptic  flat  plate  at  kS°  yaw  is  studied. 

Unlike  the  consistently  good  agreement  shown  in  the  preceding  two  figures, 

Eq.  (7-2)  in  the  present  theory  yields  a span  load  (in  thin  full  line) 
appreciably  lower  than  that  by  the  panel  method  (in  open  circles)  on  the 
forward  (left)  wing  panel.  It  must  be  pointed  out,  however,  that  even 
in  this  instance,  Eq.  (7-2)  is  still  far  better  than  the  uncorrected  strip 
theory  (in  dashes).  Its  usefulness  can  be  improved  considerably  further  by 
recasting  Eq.  (7-2)  into  a rational  fraction  similar  to  Prandtl's  original 
form,  Eq.  (7.3).  Namely: 

(S.L)  = (S.  L.)m  /[2-(S.L.)/(S.L)m]  (7.5) 

where  the  subscript  od  refers  to  4?,  -*■  oo  , noting  that  [(S-L.)  /(S.L.)^-  ij 
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belongs  to  O(sC)*  The  result  of  this  conversion  (in  heavy  full  line) 
is  shown  in  Figure  6 where  the  improvement  in  the  agreement  with  the  panel 
method  is  seen  to  be,  indeed,  significant. 


An  alternative  procedure  is  to  transform  the  wing  upwash,  say  V„  * A,  Vi  , 
instead  of  the(S.L),  l nto  V./[, V./ V.  3 ■ 1 ts  merit  remainS  t0  be  Stud'ed 
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gure  1 - Coordinates  and  Notations 
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Figure  k - Span  loading  of  an  elliptic  flat  plate  with  a ratio  of  unyawed 

span  to  root-chord  ^3  =16.78  and  a straight  b0%  chord  line  at 
Lc o • 


Figure  5 - Span  loading  of  an  ESP  flat  plate  with  a ratio  of  unyawed  span 
to  root-chord  = 16.78  and  a straight  40-6  chord  line  at  45 


Figure  6 - Span  loading  of  an  elliptic  flat  plate  with  a ratio  of  unyawei 
span  to  root-chord  =5  and  a straight  50%  chord  line  at  30° 
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